Modeling the Fermi arc in underdoped cuprates 
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Angle resolved photoemission data in the pseudogap phase of underdoped cuprates have revealed 
the presence of a truncated Fermi surface consisting of Fermi arcs. We compare a number of proposed 
models for the arcs, and find that the one that best models the data is a d-wave energy gap with a 
lifetime broadening whose temperature dependence is suggestive of fluctuating pairs. 
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I. INTRODUCTION 

It is well established that cuprates possess a supercon- 
ducting phase with an order parameter that has d-wave 
symmetry^ and for hole-doped materials, this phase ex- 
ists over a range of doping above 5%. It is also well 
established that at very low dopings, the material is an 
antiferromagetic Mott insulator. Connecting these two 
states is an unusual phase known as the pseudogap, the 
nature of which is still being debated. 2 It is felt by many 
that the proper identification of this phase will have a 
major impact on the ultimate 'mechanism' for pairing in 
cuprate superconductors. 

Angle resolved photoemission spectroscopy (ARPES) 
reveals the presence of a truncated Fermi surface in the 
pseudogap phase^^ In a study of the pseudogap versus 
temperature^ this truncated Fermi surface was denoted 
as a 'Fermi arc'. The arc was shown to be intermediate 
between the d-wave node of the superconductor and the 
complete Fermi surface of the normal state. Moreover, 
the arc appears to form by a closing of the energy gap 
of the superconducting state as the temperature is raised 
above T c . Off the arc, in the 'pseudogapped' (antin- 
odal) region of the Brillouin zone, the spectral gap ap- 
pears instead to fill in with temperature. This filling in 
effect is also seen in c-axis conductivity data£ and is con- 
sistent with the thermal evolution of the specific heat. 8 
This 'dual' nature of the energy gap is suggestive of a 
'two gap' scenario where a 'superconducting' gap resides 
on the arc and a 'pseudogap' resides off the arc. Such 
a two gap picture was proposed by Deutscher^ and re- 
cent Raman^ ARPES and STMH data have been 
offered in its support. On the other hand, even for un- 
derdoped samples, the gap function below T c seems to 
be more or less d-wave likeJ^ This conundrum of having 
a single gap below T c transforming into a dual gap above 
T c was stressed sometime ago.— 

Recently, a very detailed temperature and doping 
study of the energy gap above T c was done by Kanigel 
et al*^ They found that the length of the arc scales as 
T/T*, where T*, the temperature at which the spectral 
gap 'fills up' in the antinodal region of the zone, strongly 
increases with underdoping. 17 As a consequence, the an- 



gular anisotropy of the pseudogap looks more and more 
like a d-wave gap as the temperature is lowered relative 
to T*. This finding is supported by thermal conductiv- 
ity data, which indicates that the d-wave dispersion of 
the superconducting state at low temperatures survives 
when the doping is reduced into the pseudogap stated It 
is also consistent with recent ARPES and STM data on 
the stripe ordered phase of La7/sBa 1 / 8 Cu04, which indi- 
cates a d-wave like gap anisotropy at low temperatures 
(but above T c )™2, More recently, the study of Kanigel 
et al. has been extended to below T cf 2£ where it was 
found that the arc collapses to a node within the resis- 
tive width of the transition, with a simple d-wave like 
gap below T c . These recent studies bring into question 
the 'two gap' picture. 

Primarily motivated by the ARPES data, a wide range 
of models have been proposed to explain the Fermi arc. 
Basically, these models can be grouped into two cate- 
gories. In the first, the pseudogap is associated with 
a q=0 instability. Most of the models in this category 
have the pseudogap as a precursor to the superconducting 
gap, and involve pair formation with the absence of long 
range phase order . 21 : 22 ' 23 These models have been ex- 
tended to describe the arc by explicitly invoking vortex- 
like excitation a 24 : 25 : 26 : 27 as revealed by measurements of 
the Nernst effect^ 8 - The node of the d-wave dispersion 
is broadened into an arc by a combination of lifetime 
broadening as well as Doppler shifts of the single par- 
ticle states due to the vortices. There are, though, q=0 
theories which do not involve superconductivity. One ex- 
ample is the model of Varma and Zhu,-22 which involves 
circulating currents within a Cu02 plaquette (and thus 
has the same periodicity as the unit cell) . In this case, the 
gap function has a 'd 2 ' anisotropy. Another example is 
the 'nodal nematic' phase of Kim et al., 30 where the node 
is displaced by a nematic order parameter rather than a 
vortex Doppler shift. The final example we mention is 
the model of Wen and Lee 3 -^ where the node is displaced 
in energy due to staggered flux phase correlations. In 
fact, a rich variety of behavior has been predicted within 
the general context of resonating valence bond (RVB) 
theories.- 3 ^ 3 ^!^ 

The second category involves a non-zero q vector. This 
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category ranges from models based on a precursor spin 
density wave ) 36 i 37 charge density wave* 3 ^ stripes^ flux 
phases,— or orbital currents.— In the case of fluctuating 
order ) 32 i 41 ' 42 i 43 the non-zero q vector is not as obvious 
in the excitation spectrum. Those scenarios involving a 
(tt, tt) wavevector possess small hole pockets centered at 
(7r/2,7r/2) where the intensity is reduced on one side of 
the pocket due to the amplitude factors which mix the 
states differing by q. Related models are those where the 
Luttinger surface (surface of zeros of the single particle 
Greens function) differs from the Fermi surface i 44 i 45 ' 46 i 47 
In this case, the Fermi surface is truncated where it 
crosses the Luttinger surface. In a more general 2kp 
context, the flat parts of the Fermi surface which reside 
in the antinodal region of the zone can be eliminated by 
nesting ; i leaving a residual arc. 

In this paper, some of these scenarios will be addressed 
in the context of the ARPES data. In Section II, several 
non-zero q scenarios, where for simplicity long range or- 
der is assumed, will be analyzed. These scenarios typi- 
cally lead to (a) Fermi arcs whose length is T indepen- 
dent, (b) deviations of the arcs from the underlying Fermi 
surface, (c) energy gaps which are not centered symmet- 
rically about the Fermi energy, and (d) shadow bands. 
We argue that there is no evidence for these effects in 
ARPES and tunneling data, at least in the mildly under- 
doped region. In Section III, we turn to the q=0 solu- 
tions. We find that the scenario most consistent with the 
data is one where the node remains along the zone diago- 
nal and at the Fermi energy. The temperature evolution 
of the arc above T c is consistent with lifetime broadening 
of the node, though the data also indicate a distortion of 
the d-wave gap anisotropy with temperature. In Section 
IV, we offer some conclusions, and suggest future ARPES 
experiments that could further differentiate between the 
various models for the Fermi arc. 



II. NON ZERO Q SCENARIOS 
A. Commensurate density wave 

These scenarios assume a q vector of (tt, tt) with 
an energy gap that is either isotropic, or has d-wave 
symmetry— The secular matrix is of 2 by 2 form, and the 
Greens function associated with the wavevector k in the 
presence of simple elastic broadening, T, can be written 
as: 
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FIG. 1: (Color online) (a) Spectral intensity at zero energy 
versus k x ,k y , and (b) versus energy for k x = 1 for several k y , 
for the d-density wave 

models The black curve in (a) is the 
normal state Fermi surface. Zone dimensions for all figures are 
in ir/a units, and energies are in meV. For all figures (unless 
otherwise noted), Ao=50 meV and T=25 meV. 



We have looked at several cases, with various dispersions, 
Ck , including some with bilayer splitting, and several dif- 
ferent forms for . For brevity, we present results using 
for tk the tight binding dispersion of Norman et awS and 
a d-density wave garj 4 ^ = -^(cos(fc x ) — cos(fc y )). 

In Fig. la, we present the intensity plot of the spectral 
function (imaginary part of Gk) in the 2D zone for u) = 0. 
At the simplest level, one indeed finds an arc. But there 
are several details worth pointing out. First, the ends of 
the arc turn away from the underlying Fermi surface of 
the normal state. This is due to the fact that the zero 
energy contour traces out a pocket centered at (tt/2,tt/2), 
the back side of which is suppressed by the coherence 
factors (the prefactor of each term in Eq. I). Second, 
there is a strong suppression of the intensity at the 'hot 
spots' - the points where the normal state Fermi surface 
(ek = 0) crosses its (tt, tt) displaced image. This can 
be related to the 'Luttinger surface' effect mentioned in 
the Introduction. To see this, we note that the Greens 
function in this model can be rewritten as: 
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The 'gap' self-energy (the last term of this equation) 
diverges when u> = tk+q in the absence of broadening 
(r = 0). Thus the (tt, tt) translated image of the nor- 
mal state Fermi surface (ek+q — 0) is the Luttinger sur- 
face, and therefore the zero energy intensity is suppressed 
when the normal state Fermi surface crosses this surface. 
Finally, there is weaker intensity centered around the 
(tt, 0) points which will be suppressed as Ao increases. 
To investigate this further, in Fig. lb, we show the spec- 
tral function for several k points along the (tt, 0) — (tt, tt) 
direction. One clearly sees that the spectral function 
has a minimum value that sits at negative energy. At 
k = (tt, 0), it is obvious from Eq. 2 that this minimum 
value occurs at u) = e^, which is -34 meV for this disper- 
sion. This asymmetry in energy is obviously enhanced 
for dispersions where e^o is deeper in energy. 
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FIG. 3: (Color online) (a) Spectral intensity at zero energy 
versus k x ,k y , and (b) versus energy and k y for k x = 1, for the 
model of Yang et al£^ 



ident. This is consistent with tunneling data as well 51 
where the minimum in the tunneling conductance is at 
zero bias, even in the pseudogap phase. 



FIG. 2: (Color online) (a) Experimental energy distribution 
curves (EDCs) for optimal doped Bi 2 Sr 2 CaCu208 (Bi2212) 
around the underlying Fermi surface in the pseudogap phase 
(T=140K) divided by a resolution broadened Fermi function. 
The bottom curve is at the node and the top curve at the 
antinode. The data set is the same as in Fig. lb of Ref.lla (b) 
Zero energy intensity from (a) as a function of the Fermi sur- 
face angle, (j>, where 4> — 0° corresponds to the antinode and 
= 45° to the node, (c) Zero energy intensity versus k x ,k y 



(the data were reflected relative to k x 
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For (b), the 



intensities were obtained by subtracting a background EDC 
(obtained from an unoccupied k), then normalizing this sub- 
tracted intensity by its energy integrated weight. This was 
designed to minimize the effect of the photoemission matrix 
elements. This was not done in (c) in order to demonstrate 
that the raw data show no indication for any deviation of the 
arc from the underlying Fermi surface (black curve). 



In relation to the experimental data, we note the fol- 
lowing issues with this model, which are generic to mod- 
els based on a finite q order parameter. First, there is 
no natural way to generate an arc whose length is pro- 
portional to the temperature. Second, there is no ev- 
idence from ARPES for a 'turn in' of the ends of the 
arc away from the underlying normal state Fermi surface 
(Fig. 2c). Third, ARPES is consistent with spectral func- 
tions which either have a maximum (arc) or minimum 
(antinodal region) at zero energy along the underlying 
Fermi surface. We demonstrate this in Fig. 2a, where 
data in the pseudogap phase along the underlying Fermi 
surface is plotted. These data are the same as used to 
construct Fig. lb of Ref. [TH, but instead of 'symmetriz- 
ing' the raw data as was done there (which implicitly 
assumes a maximum or minimum at zero energy) , we di- 
vide the data by a resolution broadened Fermi function. 
The clear maxima at zero energy along the arc, and the 
minima at zero energy away from the arc, are quite ev- 



B. Differing Luttinger surface 

These scenario o 44 ' 45 ' 46 ' 47 are related to the ones just 
discussed. For discussion purposes, we look at the re- 
cently proposed model of Yang, Rice, and Zhang. 44 In 
this case, the Greens function is 
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where is just the near neighbor term of the tight 
binding dispersion (A& has the same form as the d 
density- wave case) . Note that if ek only had a near neigh- 
bor term, then at half filling, this model would be equiva- 
lent to the d density- wave model. The similarity of these 
two models can be seen in Fig. 3, where we show the zero 
energy intensity plot in the 2D zone, as well as the inten- 
sity versus u> for k along (tt, 0) — (7r,7r). Again, note the 
pronounced suppression of the intensity at the 'hot spots' 
in Fig. 3a, which is not evident in the data (a plot of 
the experimental zero energy intensity around the Fermi 
energy is shown in Fig. 2b), as well as the pronounced 
asymmetry of the energy gap relative to the Fermi energy 
in Fig. 3b. And, as with the d density-wave model, there 
is no obvious mechanism to obtain an arc proportional 
to T. 



C. Nesting density wave 

These scenarios assume a q vector which nests the 
antinodal points of the 2D Fermi surface. Two approx- 
imations were analyzed. In the first, a single q vec- 
tor along q yi q = (0, — q), was used in the first octant 
(bounded by (0, 0) — (tt, 0) — (tt, tt) — (0, 0)) of the square 
lattice zone (a 2 by 2 secular equation), the result of 
which was then reflected to the other octant. The orien- 
tation of q was designed so as to connect the antinode at 



4 




100 



LU 



100- 




-100 



0.0 



ky 



0.5 



0.0 



ky 



0.5 



FIG. 4: (Color online) (a) Spectral intensity at zero energy 
versus k x , k y for (a) the 2 by 2 secular equation and (b) the 3 
by 3 secular equation approximations, for the antinodal nest- 
ing model with q = (0.36, 0). 



FIG. 5: (Color online) (a) Spectral intensity versus energy 
and k y near the arc tip (k x — 0.6) for (a) the 2 by 2 secular 
equation and (b) the 3 by 3 secular equation approximations, 
for the antinodal nesting model with q — (0.36, 0). 



(ir, q/2) with the one at (w, —q/2). The equation for the 
Greens function is the same as in Eqs. 1 and 2, except 
that Afc in this case was taken to be isotropic. 

In the second approximation, a 3 by 3 secular equa- 
tion is separately solved for q vectors oriented respec- 
tively along q x , q = (q, 0) and q = (—q, 0), and along q y , 
q = (0,g) and q — (0, — q), in the first quadrant of the 
zone, and then the two results are averaged (represent- 
ing averaging over two domains). The unaveraged G\. is 
given by 
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where by i + 1 and i + 2 we mean modulo 3. The Ei are 
given by solving the appropriate cubic equation and can 
be written as 
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In Fig. 4, we show the zero energy intensity in the 2D 
zone for the two approximations. Again, a clear arc is 
seen, with extra structure that can be attributed to the 
reduced intensity (due again to the coherence factors) 
of the 'shadow' bands. This is particularly true in 4b 
where more shadow bands occur. A similar situation 
would occur if one had 'checkerboard' order (this would 
be obtained by solving a 5 by 5 secular matrix associated 
with a 'double q' structure). 

A significant difference from the previous cases is the 
origin of the arc. In the previous cases, the arc is due to 
the Fermi energy cutting across the lower of the two en- 
ergy bands. In essence, the energy gap is centered above 



the Fermi energy for k vectors from the node to the 'hot 
spots', and it is centered below the Fermi energy for k 
vectors from the 'hot spots' to the antinode. But in this 
'antinodal' nesting case, it is the reverse situation. In 
the 2 by 2 approximation, the arc is formed from the 
Fermi energy cutting across the upper of the two bands. 
This is particularly evident near the arc tip, as shown in 
Fig. 5. In essence, the energy gap is centered below the 
Fermi energy for k along the arc. For momenta near the 
arc tip, one would find a minimum in the spectral func- 
tion at a negative energy (as in Fig. lb). This effect is 
not evident, though, in the ARPES data. And again, as 
the arc tip is associated with where the underlying nor- 
mal state Fermi surface intersects the density wave zone 
boundary (in the first octant, this would correspond to 
k y = q/2), there is no natural mechanism for an arc pro- 
portional to temperature. As discussed by McElroy, 49 
this would require a 'two gap' scenario, where the den- 
sity wave gap would wipe out the antinodal parts of the 
Fermi surface, and then a second gap would wipe out the 
remaining arc with reducing temperature. Despite the 
attractiveness of such scenarios in regards to some ex- 
perimental data^ ' 10 i n ' 12 i 13 a definitive signature of this 
density wave gap would be to observe the shadow bands 
evident in Fig. 4 and the asymmetry of the gap relative 
to the Fermi energy evident in Fig. 5. So far, we have 
found no evidence for either of these effects^ 2 - 



III. ZERO Q SCENARIOS 

A. Energy displaced node 

The RVB model of Wen and Lee 3 -! is based on incor- 
porating both the effect of a d-wave gap in the particle- 
particle channel and a staggered flux phase gap in the 
particle-hole channel. An ansatz for the Greens function 
in this model that makes it of the same form as the earlier 
cases we studied is 
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FIG. 6: (Color online) (a) Spectral intensity at (a) zero energy 
and (b) -25 meV versus k x , k y for the model of Wen and Lee 
with /i sh =-50meV .— 

The effect of fj, s h is to move the d-wave node off the Fermi 
energy, and the sign of fx s h is chosen to be negative so that 
the node is above the Fermi energy. The result is an arc 
at zero energy that is tied to the underlying Fermi surface 
(Fig. 6a). This is a positive feature of this model. If \fi s h\ 
were proportional to T for some (unknown) reason, this 
could also account for the temperature evolution of the 
arc. The major problem with this model, though, is that 
as one goes to positive energy, the arc would shrink in 
size, and as one goes to negative energy, the arc would 
expand in size. This effect is not evident in the ARPES 
data in the pseudogap phase, though, which seems to be 
more or less consistent with arcs which are independent 
of energy up to energy scales of order Ao .— And in this 
model, the energy gap is always centered at an energy 
above the Fermi energy for all momentum cuts in the 
zone, an effect not visible in Fig. 2a. 

B. d-wave pairs plus lifetime broadening 

The simplest model in this class is equivalent to the 
one just described with = 

G^=cj-e k + iT- — ^— , (9) 
lu + e k + iT 

The spectral function for finite T traces out an 'arc', as 
shown in Fig. 7a. And the energy gap is centered at the 
Fermi energy, as shown in Fig. 7b. That is, the gap is 
tied to the Fermi energy and the Fermi surface, consistent 
with experiment. In Fig. 8a, we plot the evolution of 
the spectral function on the Fermi surface (e^ = 0) for 
this model, and in Fig. 8b the angular anisotropy of the 
spectral gap (half the peak to peak separation) . 

Gapped and ungapped spectra on the Fermi surface 
(Fig. 8a) are obviously controlled by the sign of the sec- 
ond derivative of the spectral function with respect to uj 
at uj — 0. The condition that this second derivative is 
zero is r = i/3A k ■ Assuming a simple d-wave gap of the 
form Afe = Ao cos(2</>) where 4> is the Fermi surface angle 
measured relative to the antinode, one then obtains for 
the position of the arc tip (j>o = 0.5 cos _1 (r/\/3Ao). T* 
would then be the condition that T(T) = V3A (T). 




0.0 1.0 0.0 0.5 

kx ky 

FIG. 7: (Color online) (a) Spectral intensity at zero energy 
versus k x ,k y , and (b) versus energy and k y for k x = 0.6, for 
the d-wave pair model. 
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FIG. 8: (Color online) (a) Spectral intensity around the Fermi 
surface for the d-wave pair model. The top curve is at the 
node, the bottom curve at the antinode. (b) Spectral gap 
(half the peak to peak separation) versus the Fermi surface 
angle. The dashed curve corresponds to T=0. 

In Fig. 9a, we show the variation of the arc length 
with r. This variation is consistent with experiment, as 
shown in Fig. 9a, if one assumes that r cx T and An is 
a constant in T (similar plots are shown in Refs. 0JH3). 
This linear variation of the arc length with T is a natural 
consequence of the linear variation of A k with (f> around 
the node. The upturn of the arc length at larger T is 
due to the quadratic dependence of the energy gap with 
4> about the antinode. 

The theory of Varma and Zhu 29 is similar except that 
Afc is taken to be the square of the d-wave gap. As they 
point out, this fits the angular anisotropy of the param- 
eter Afe in the pseudogap phased better than the simple 
d-wave model, as can be seen in Fig. 9b. On the other 
hand, the arc length variation with T is more consis- 
tent with the simple d-wave form, as shown in Fig. 9a, 
though we remark that Varma and Zhu were able to ob- 
tain a much better fit to the arc length by allowing a 
self-energy with a more sophisticated frequency and tem- 
perature dependence^ 



IV. CONCLUSIONS 

In regards to the 'non zero q' scenarios, there are sev- 
eral ways that experiment could address this question. 
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FIG. 9: (Color online) (a) Arc length versus V for a d-wave 
gap and a d 2 gap (data from Ref. |T(i '■ . (b) Experimental 
spectral intensity at zero energy versus Fermi surface angle 
(as in Fig. 2b) versus fits that assume either a d-wave gap 
(cos(2<^))) or a d 2 gap (cos 2 (20)). 



Definitive evidence would be finding a departure of the 
arc from the underlying normal state Fermi surface, ev- 
idence for an energy gap which is asymmetric in energy 
relative to the chemical potential, or the existence of 
shadow bands (i.e., bands displaced from the main band 
by the wavevector q). Other evidence would be the ex- 
istence of intensity suppression at 'hot spots' (where the 
Fermi surface would cross the Luttinger surface), as has 
been observed by ARPES in electron doped cuprates^ 
In regards to the 'zero q' scenarios, the simplest the- 
ory consistent with the data appears to be a d-wave gap 
with an inverse lifetime that is proportional to T. There 
are, though, some limitations of this model. The data of 
Ref. [l6| were actually fit with a form more general than 
that of Eq. 9^ 
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This 'two lifetime' model has the advantage of being 
able to describe a broad spectral function (Ti) but with 
a sharp leading edge gap (r ) as indicated by ARPES 
data in the pseudogap phase4&££ It has since been ex- 
tended to include a more general frequency dependence 
for the self-energy.— The presence of two lifetimes may 
seem unusual, as this does not occur, for instance, in 
the standard Eliashberg treatment of strong coupling 
superconductors^ The motivation in Ref. |56| was that 
Ti denotes the interaction in the particle-hole channel, 
whereas Tq denotes that in the particle-particle channel. 
If the dispersion of the single particle states is ignored, 



a calculation of fermions interacting with pair fluctua- 
tions leads to Tq proportional to T — T c £^ In contrast, in 
two dimensions, a dependence proportional to \JT — T c 
is obtained*^ Consideration of vortex excitations 2 ^ leads 
instead to linear T behavior at high T . In general, one 
would expect linear T behavior at high temperatures, 
since this occurs for any model of fermions interacting 
with bosons^! - but with a collapse to zero at T c since the 
inverse pair lifetime should vanish in the ordered state. 
This general T dependence not only naturally describes 
the linear T variation of the arc length^ it can also 
account for the 'filling up' of the gap in the antinodal 
region. 5 - It also explains why the arc length collapses to 
zero within the resistive width of the transition. 20 Re- 
gardless, our experience has been that two lifetimes are 
necessary to properly model the data i 16 ' 56 i 57 A complete 
description of modeling based on Eq. 10 in regards to 
the ARPES data is beyond the scope of the present pa- 
per and will be left for a future study. 

The fits presented in Ref. |l6| also indicated that the 
gap anisotropy changed with temperature, and that this 
effect could not be described by lifetime broadening of 
the zero temperature gap, though it should be remarked 
that the actual value of is difficult to extract once 
the broadening significantly exceeds A/.. Whether such 
changes in anisotropy (and in particular, a region around 
the node where A& is identically zero as indicated by 
these fits) can be described by pair breaking within a 
strong-coupling Eliashberg context remains to be seen. 
Certainly, careful measurements of the gap anisotropy at 
different temperatures and dopings would help to better 
differentiate models based on a d-wave gap from more 
novel ones, such as that of Varma and Zhu^ that have 
a fundamentally different gap anisotropy. 

In conclusion, we believe that a model of a d-wave gap 
with a temperature dependent lifetime consistent with 
fluctuating pairs gives the simplest description of ARPES 
data in the pseudogap phase. Future experimental work 
should be aimed at further differentiating between vari- 
ous proposed models for the Fermi arc, as well as using 
this information to address other data, such as transport. 
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